Semi-relativistic meson exchange currents in (e, e') and 

(e, e'p) reactions 



J.E. Amaro 1 , M.B. Barbaro 2 ' 3 , J. A. Caballero 3 , F. Kazemi Tabatabaei 1 

1 Departamento de Fisica Moderna, Universidad de Granada, E-18071 Granada, SPAIN 

2 Dipartimento di Fisica Teorica, Universita di Torino and INFN, Sezione di Torino 
Via P. Giuria 1, 10125 Torino, ITALY 

3 Departamento de Fisica Atomica, Molecular y Nuclear 
Universidad de Sevilla, Apdo. 1065, E-41080 Sevilla, SPAIN 



Abstract 

Electron-induced one-nucleon knock-out observables are computed for moderate to high mo- 
mentum transfer making use of semi-relativistic expressions for the one-body and two-body meson- 
exchange current matrix elements. Emphasis is placed on the semi-relativistic form of the A-isobar 
exchange current and several prescriptions for the dynamical-equivalent form of the A-propagator 
are analyzed. To this end, the inclusive transverse response function, evaluated within the context 
of the semi-relativistic approach and using different prescriptions for the A-propagator, is com- 
pared with the fully relativistic calculation performed within the scheme of the relativistic Fermi 
gas model. It is found that the best approximation corresponds to using the traditional static A- 
propagator. These semi-relativistic approaches, which contain important aspects of relativity, are 
implemented in a distorted wave analysis of quasielastic (e, e'p) reactions. Final state interactions 
are incorporated through a phenomenological optical potential model and relativistic kinematics 
is assumed when calculating the energy of the ejected nucleon. The results indicate that meson 
exchange currents may modify substantially the TL asymmetry for high missing momentum. 
PACS: 25.30.Fj, 14.20. Gk, 24.10.Jv, 24.10.Eq, 21.60.Cs 

Keywords: Nuclear reactions, Exclusive electron scattering, Inclusive electron scattering, Meson 
exchange currents, Relativistic Fermi Gas. 



1 Introduction 



Electron induced one-nucleon emission reactions near the quasielastic peak are clearly dom- 
inated by one-body (OB) dynamics since such kinematics roughly corresponds to the virtual 
photon interacting directly with a bound nucleon [1, 2, 3]. Under these conditions it has 
been possible to analyze a considerable amount of data for a variety of nuclei and to study 
many aspects of the reaction mechanism. Most of these studies have been based on the stan- 
dard distorted- wave impulse approximation (DWIA), where the OB current matrix elements 
are computed using single-particle wave functions obtained as solutions of the Schrodinger 
equation with phenomenological potentials. Recently [4, 5, 6, 7], efforts have been made to 
develop a fully relativistic description of the process in order to describe consistently the 
high momentum and energy transfer region [8]. This constitutes the basis of the relativistic 
distorted- wave impulse approximation (RDWIA). 

Within the scheme of non-relativistic approaches, the role played by correlations beyond 
the mean field and their impact on the overlap functions and spectroscopic factors has 
been investigated in several works. Many-body calculations [9, 10] including short-range 
correlations (SRC) of central type only, have shown a small effect over the quasi-hole overlap 
functions and, accordingly, over the (e, e'p) cross section. In [11] the SRC were shown not 
to modify substantially the mean field (e, e'p) results at high momentum and low excitation 
energy. In recent works a significant effort has been made to improve the analysis by including 
the tensor and spin-isospin channels [12, 13] as well as long-range correlations [14, 15]. It is 
important to point out that the extraction of spectroscopic factors from the analysis of (e, e'p) 
experiments is still not free from ambiguities. This would require an accurate knowledge of 
the reaction mechanism. In this sense, a RDWIA calculation, compared to DWIA, gives rise 
to a different quenching of the (e, e'p) cross section, hence leading to different spectroscopic 
factors [4] (larger for RDWIA). 

The two-body meson exchange currents (MEC) may also play a significant role in the 
description of electron scattering observables since they are connected to the nuclear corre- 
lations by the continuity equation. In fact the interplay between correlations and MEC in 
(e, e'p) reactions is far from trivial; rigorously the concept of overlap function is not enough 
to describe the process in presence of two-body current operators, and therefore the effect 
of correlations in presence of MEC may be not simply parameterizable by a spectroscopic 
factor. This idea is supported by calculations of the inclusive transverse response for nu- 
clear matter within the correlated basis function perturbation theory [16], that have shown 
a significant effect due to MEC, contrary to the much smaller effect usually found for uncor- 
rected calculations [17]. Before a complete calculation of the exclusive response functions 
including MEC within a sophisticated correlated model of the reaction be attempted, it is 
necessary to calibrate different approaches to the calculation in uncorrelated models, as a 
first step beyond the DWIA. 

Thus in this work we restrict our attention to the mean-field (uncorrelated) model and 
evaluate the role played by the MEC in quasielastic inclusive and exclusive electron scat- 
tering reactions. In the case of exclusive (e, e'p) processes, a theoretical evaluation of MEC, 
including the usual pion-in-flight (P), contact (C) and Delta-Isobar (A), has been presented 
in refs. [18, 19, 20]. In [20] substantial differences were found with respect to the results 
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of ref. [19], particularly for the interference TL response. The analysis in [19], which treats 
final state interactions (FSI) through a real potential, was extended to higher momentum 
transfer in refs. [21, 22] where some relativistic corrections were also included in the OB 
current operator. Instead the standard non- relativistic MEC operators [23, 24] were used 
except for a modified "dynamical" A-propagator that depends on the invariant energy y/s 
of the A. More recently, the theoretical model of [18] has been refined in [25] modifying 
also the A-propagator in order to include the A-invariant mass. The use of a proper "dy- 
namical" propagator in the A-current has been discussed at length, concerning mainly the 
case of two-particle emission reactions [26, 27, 28]. One of the central questions is how to 
determine y/s in a way that could be easily implemented in existing models for (e, e') and 
(e, e'p) reactions. 

In general the appropriate choice of y'i will depend on the specific dynamical model for 
the reaction mechanism. Hence a proper description of the reaction dynamics through an 
adequate choice of the currents and kinematics is necessary. For high momentum transfer 
relativistic effects play a significant role and therefore relativity should enter in the descrip- 
tion of kinematics and/or current matrix elements. Within a fully relativistic model, MEC 
effects in (e, e'p) have been computed recently in [29], where only the contact (or "seagull") 
current has been considered. Due to the complexity of the fully relativistic two-body cur- 
rents, particularly the A one, in this work we start by developing reasonable semi-relativistic 
expressions for all the two-body currents (P, C and A) that can be easily implemented in 
existing non-relativistic descriptions of the reaction mechanism. These semi-relativistic cur- 
rents retain important aspects of relativity and hence they may be used to describe properly 
recent experiments performed at high momentum and energy transfers. 

In this paper we discuss the one-particle emission sector, introducing a recently developed 
semi-relativistic (SR) approximation for the three MEC operators which, in conjunction 
with the SR form of the OB current derived in [30] (see [31] for a recent review on these 
expansions), makes it possible to evaluate quasielastic (e, e') and (e, e'p) observables. The SR- 
OB current was applied in inclusive (e, e') [32], and for coincidence (e, e'p) [6, 10, 33, 34, 35] 
reactions. The SR expressions for the MEC operators were obtained in [36] in the case of 
the pionic and contact currents. In [31, 37] a comparison between the SR-MEC predictions 
for the inclusive transverse response and the exact RFG results is presented. 

One of the goals of this work is to derive a semi-relativistic expression for the A-exchange 
current and compare its prediction for the inclusive transverse response function with the 
exact relativistic result obtained within the RFG model of ref. [38] . This comparison, carried 
out in the inclusive channel (section 2), makes it possible to test the reliability of the different 
prescriptions of the "dynamical" A-propagator to be used in exclusive (e, e'p) processes 
(section 3). Once this is settled, the impact of the relativistic MEC over the exclusive 
response functions is evaluated for moderate to high momentum transfer. Finally, in section 4 
we draw our conclusions. 
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Figure 1: Feynman diagrams contributing to the two-body current with one pion-exchange. 
Contact (C) (a,b), pionic (P) (c), and isobar (A) (d)-(g) are included in this work. 



2 Inclusive transverse response 

The impact of meson exchange currents on inclusive electron scattering has been tradi- 
tionally investigated in a non-relativistic framework and relativistic corrections have been 
implemented making use of a power expansion in the momenta of the nucleons and of the 
exchanged photon [39, 40, 41]. In refs. [31, 38] a semi-relativistic analysis of the MEC (P, 
S and A) currents, represented in fig. 1, has been performed in the context of the RFG 
with the aim of finding simple prescriptions for implementing relativistic effects into non- 
relativistic calculations. Within the RFG model, it is possible to perform a fully relativistic 
calculation of the inclusive response functions, using exact relativistic kinematics, currents 
and propagators. Hence the inclusive transverse response function of nuclear matter can be 
used to test the quality of the semi-relativistic MEC, to be applied in the next section to 
semi-inclusive processes. 

In this section we focus on the A current, represented by diagrams (d-g) of fig. 1, putting 
emphasis on the "dynamical" treatment of the A propagator, an aspect of the problem which 
was neglected in the SR model of [38]. Restricting our attention to the case of one-particle 
emission processes, the contribution of the A current to the transverse nuclear response 
function is obtained as the interference between the A and the OB currents. In the Fermi 
gas model it reads 



R 



a-OB 



3Z 



2Re / d 3 h5(E p -E h - cu)Af ph Ti \)° B (p, h)* • j£(p, h) 



h<kp 



3 



where p = h + q is the momentum of the ejected particle, is the transverse (perpendicular 



to q) component of j and kp is the Fermi momentum. In the RFG model E^ = y k 2 + m N , 
whereas in the non-relativistic (NR) Fermi gas, E^ = m^ + e^ = + k 2 /(2mjv). Moreover, 
the factor A/" p h, arising from the spinor normalization, is m 2 N j '(E p E h ) for RFG and 1 in a 
NR model. Finally j^(p,h) (which implicitly include spin and isospin indices) are the OB 
or A particle-hole currents [31]. 

The non relativistic two-body A current used in this paper has been obtained from the 
relativistic one in the limit of small momenta by following the procedure sketched in ref . [38] ; 
it can be written in momentum space as 

jnr(p'l ; P2 ) Pl>P2) 

= {Ga(Pi + Q)q x [-k 2 x <r« + 2^k 2 \ [2rf - i[r^ x r (2) ]J 

i (2) 

+ G A (P[ - Q)q x [k 2 x tr<» + 2ik 2 l [2rf + i[r<» x r (2) ]J } ff-^ 

k 2 + 

+ (1 <— > 2) , (2) 

where Pi, P[ are the four- momenta of the initial and final nucleons defined in fig. 1, k; = 
p'j — pi is the momentum transferred to the i-th nucleon and Ga(P) is the non-relativistic 
(dynamical) A propagator, which will be later discussed in detail. We use the following 
values for the coupling constants: Gi = 4.2, / 2 /47r = 0.079, f w NA = 2.24. 

In [38] a more general relativistic current was considered, containing additional couplings 
G 2 , G 3 and dependence on the so-called off-shell parameters zi, z%, z 3 . Therein it was shown 
that the terms G 2 and G 3 have a weak impact on the transverse response for q below 1 GeV 
and that this response displays a small dependence upon the off-shell parameters compared 
with its sensitivity to y/s, which will be discussed at the end of this Section. Therefore in this 
work we use the Peccei lagrangian, which neglects G 2 and G3 and corresponds to Zi = —1/4. 

The MEC ph matrix elements are obtained by summing over all the single-particle states 
occupied in the ground-state Slater determinant [31, 38]. In the case of nuclear matter only 
the exchange term, represented by the many-body diagrams of fig. 2, survives, due to spin 
and isospin symmetries. Moreover, it can be proved that the exchange A terms represented 
by diagrams (f) and (g) are also zero in nuclear matter. In the case of finite nuclei there 
exists however a contribution coming from the direct A term, which, although small, has 
been included in the calculations presented in next section. 

Using eq. (2) the particle-hole matrix element of the non-relativistic A current in nuclear 
matter results 

4 Gi f wNA f r d 3 k 



•A 1 i\ -x 4 Gl J-* NA J f d k an i\ 



x {G A (k + q)B(h - k) SpSh + G A (k - q)B(p - k) SpSh } , (3) 
where = 1 if h is a proton and —1 if neutron, and we have defined the function 
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Figure 2: Exchange MEC diagrams contributing to the one-particle emission responses com- 
puted in this work. In the vertex-type isobar diagrams (d,e) different prescriptions for the 
dynamical A propagator are discussed in the text. 



If the k-dependence of the A-propagator is neglected, then the integral in eq. (3) can be 
solved analytically 1 (see ref. [17] for its explicit form). 

The relativistic expression for the currents can be found in [38] and yields results which 
differ significantly from the non-relativistic ones even if the momentum transfer is as low as 
500 MeV/c. This difference is due in part to the relativistic kinematics (RK), which can 
be easily implemented in the NR model by the replacement u — > u>(l + oo/2m N ) [32], with 
the exception of the electromagnetic form factors, that must be computed for the unshifted 
value of uj. Moreover, the NR results can be brought closer to the relativistic ones if the SR 
form of the current operators 

jr,5fl(P' h ) = -7T^jr,nr(P»h) (5) 

Vl + T 

is employed for a =OB,C,P,A, with r = \Q 2 \/4m 2 N . 

Eq. (5) has been obtained (see refs. [31, 36, 38] for details) by a direct Pauli reduction, 
expanding in powers of h/m N to first order the OB current [30] and to leading order the MEC 
of pionic and contact types [36, 37]. In this paper we use a similar 'factorized' expression 
for the SR A-current that was already proposed in [38]. The r-dependent factor in eq. (5) 
arises from the spinology and produces a reduction of the responses. 

1 Note that in this section we are not including 7rN form factors in order to obtain an analytical result for 
numerical convenience. In the next section they will be properly included. 
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Let us now discuss the A propagator G A , appearing in the A current, eqs. (2,3). In 
a fully relativistic theory, G A is given by the Rarita-Schwinger tensor [38]. The possible 
decay of the isobar into a Nit state is accounted for by the substitution of the A mass 
m A — > m A — |r(P 2 ) in the denominator of the A propagator, where the function r(P 2 ) is 
the width of the resonance (see [42, 43]). 

The non-relativistic version of this propagator is defined by the positive energy sector for 
spatial indices and small momenta [38] 

G U P ) - Pz^T (Sa + p7i) — G A (P) - l -e l3k a k ) . (6) 

Note that diagrams (d-g) of fig. 1 correspond to different momenta in the propagators, 
namely G A {P\ + Q) and G A {P[ — Q): these, referred to as A-excitation (G A ) and A- 
deexcitation (G A ) in [21, 25], denote the propagator for a A created after and before photo- 
absorption, respectively. 

The static limit approximation to G A consists in taking 

P + m A ^ mjv7o + m A ^ m N + m A _ 1 
P 2 — m A P 2 — m 2 A m 2 N — m 2 A m N — m A ' 

where we have used P 2 ~ m 2 N for P = Pi + Q or P = P[ — Q and q, u small, since Pi and P[ 
are the four-momenta of the initial and final nucleons. In this case, the two propagators are 
equal, i.e., G A (P\ + Q) = G A (P[ — Q) = (m N — m A ) _1 and can be factorized out in eq. (2), 
thus yielding the traditional form of the non-relativistic current [17, 23]. 

The effects introduced by relativity can be appreciated in fig. 3. Here we show the 
inclusive transverse response of 40 Ca for two values of the momentum transfer, q = 0.5 
GeV/c (top panel) and q = 1 GeV/c (bottom panel). We use Fermi momentum Icf = 237 
MeV/c and the Galster parameterization of the nucleon form factors [44]. For the A current 
we use the electric form factor of the proton. For the comparison with the non-relativistic 
result we use strong form factors F n NN — F n ^ A = 1. The fully relativistic calculation 
performed within the RFG model of ref. [38] (thick solid line) is compared with various 
non-relativistic approaches: the traditional non-relativistic model (NR), the NR model but 
including relativistic kinematics (denoted as RK), and the results corresponding to the semi- 
relativistic currents, using in addition relativistic kinematics (denoted as SRO). Apart from 
the RFG, where the Rarita-Schwinger propagator is used including the A width, the other 
results have been evaluated employing the static limit approximation for the A propagator, 
for which y/s = m N , and therefore Ta(s) = 0, since s is below the threshold, mjq + m^. Note 
the crucial role played by the relativistic kinematics and, moreover, how the SR approach, 
compared to the RK case, improves the agreement with the RFG results. 

Next we study the effect of including a "dynamical" A propagator within the SR model. 
This can be done by using different approximations for the propagators, G A (K + Q) and 
G A (K — Q), appearing in the non-relativistic current, eq. (3), under the requirement that 
the new "dynamical" propagator be independent of k: only in this case it can be easily 
introduced into non-relativistic calculations. 

The results for the A-OB transverse response computed within the SR model, and using 
different prescriptions for the A propagator, are displayed in fig. 4. As in the previous figure 
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Figure 3: Contribution of the A to the transverse response function for q = 500 and 1000 
MeV/c. The results corresponding to the relativistic Fermi gas (RFG) are compared with 
the non- relativistic (NR) Fermi gas, including relativistic kinematic (RK), and using the 
semi-relativistic (SR0) approach for the electromagnetic currents. 

and for comparison, we present the fully relativistic results (RFG) and the SR approach 
using the static limit for the A propagator (SR0). In all the cases, relativistic kinematics is 
assumed. The various approaches considered to describe the "dynamical" A propagator are 
below listed and discussed in some detail. 

1) The propagator is written, as suggested in refs. [26, 28], as 

Ga(P) = -7= \ tr n , (8) 

where s = P 2 and yfs is the available energy in the A rest system. The approximate 
values of s are obtained by neglecting the momentum k and kinetic energy ek compared 
with q and the nucleon mass m^, i.e., 

s 1 = (K + Q) 2 = (E k + Lu) 2 -(k + c l ) 2 ~(m N + Lu) 2 -q 2 (9) 
s 11 = {K-Q) 2 ~{m N -uj) 2 -q 2 . (10) 
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As we observe in fig. 4, this prescription (denoted as SRI) produces a hardening of the 
response function, that is more pronounced for high q, and under-estimates the exact 
result at the peak. 

2) This second prescription is obtained, as suggested in refs. [25, 28], by using the static 
approximation in the de-excitation diagram, fig. 2e, i.e., s 11 ~ m%, and a dynamical 
propagator in the A excitation diagram, fig. 2d. In this case we use the same value s 1 
as in eq. (9). As shown in fig. 4, this prescription (labelled as SR2) still produces a 
hardening of the response, but now the response function is larger than the relativistic 
one. We note also that the SR2 curve crosses the static SRO result precisely at the 
position of the quasielastic peak (QEP). In fact, this is determined by {tun + uj) 2 — q 2 = 
m 2 N) hence eq. (9) gives the static result, = m^. 

3) From eq. (6) we may write 

P + m A -E70 - q • 7 + rn A 1 



P 2 - m 2 A ^/s + m A - m A 



'ID 



for P = K + Q. The result given by eq. (8) is re-obtained when the first term in the 
right hand side of eq. (11) is close to one. Note that for q large, \fs can be significantly 
different from E, while one cannot neglect the term q-7. This would give an additional 
spin-dependent term and the propagator would get tangled. Instead, a possible way to 
proceed is by taking partially into account the g-dependence in the numerator in the 
way: 

n (v y rc\ m N + uj + q + m A 1 

G A {K + Q) ~ =- -= (12) 

v s + m A v s — m A 

n (is r>\ m N + uj - q + m A 1 

G A (K-Q) ~ 7 = r = , (13) 

vs 11 +m A v s 11 — m A 

where s I,n are given in eqs. (9,10). This procedure allows us to write down an ap- 
proximate diagonal expression for the A propagator, and the corresponding results are 
denoted in fig. 4 as SR3. The reduction of the response function is shown not to be as 
large as in the SRI case. 

Summarizing, from fig. 4 we may conclude that none of the above approximations pro- 
posed for the dynamical A propagator is entirely satisfactory. In fact we observe that the 
best result, compared with the fully relativistic RFG one, corresponds to the static limit 
approximation. 

To complete this discussion, we note that the role of the width Y in the SR models of 
fig. 4 is irrelevant since the approximate values taken for ^Js are far below the pole for these 
kinematics. In fact, for the prescription SRI we have from eq. (9): 



u = Js 1 + q 2 - m N . (14) 
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Figure 4: Contribution of the A to the transverse response function, using different pre- 
scriptions for the A propagator, discussed in the text 

This means that for q — 1 GeV/c and \fs^ = toa, we get u ~ 649 MeV, which is well above 
the allowed region of the QEP, while the threshold, \fs^ = + m n , is reached for uj ~ 531 
MeV, below which the A width is zero. Hence only the small tail of the A close to the 
threshold is being considered in this region of energy. 

Actually in the relativistic model things are totally different because the pole is reached 
inside the allowed energy region. Indeed the delicate energy balance in the denominator 
makes the value of the inner momentum k to play a role. Instead of (14) a better approxi- 
mation to the exact relation is represented by 

u = yV + (q + k) 2 - m^v . (15) 
In this case the pole is first reached at 

u = \f m 2 A + (q - k F ) 2 - m N , (16) 

which gives u ~ 511 MeV, close to the peak for q — 1 GeV/c. 

Thus we have considered two new prescriptions, denoted as SR4 and SR5. To help the 
reader they are shown in a separate figure (fig. 5), where they are again compared with 
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the fully relativistic calculation (RFG) and the static limit approach (SRO). These new 
prescriptions are based on the value assigned to \/s, such that the position of the pole be 
closer to the relativistic case. More precisely: 

4) One considers v^" = \/s~nn — m N, being s^n the invariant energy of the two outgoing 
nucleons, suggested in [26, 28] for two-nucleon knock-out, and apparently applied also 
to one-nucleon emission in [25]. The meaning of snn for this case is doubtful, since 
there is only one particle in the final state. However, from the A-excitation diagram, 
fig. 2d, we see that k, h are two entering momenta and p, k the two exiting ones. Hence 
the prescription SR4 in fig. 5 is based on the approximation: 

s NN ~ (P + Kf = (E h + £ k + lo) 2 - (h + k + q) 2 ~ (2m N + uf - q 2 . (17) 

In this case the pole for q — 1 GeV/c is reached at oo ~ 513 MeV, inside the QE region. 
However, from the results in fig. 5, this prescription appears to be less satisfactory than 
the previous ones. 

5) Finally, the prescription SR5 is obtained by exploiting eq. (16) that gives the right 
position at which the pole is first reached and defining 

s 1 = (m N + lo) 2 - (q - k F f (18) 
s" = (m N - lo) 2 -(q- k F ) 2 . (19) 

The results show that for q — 1 GeV/c the pole is reached at the right position close 
to the peak of the response. The change of sign of the response is a consequence of the 
change of sign of the A propagator when one crosses the pole. 

From fig. 5 we observe that it is not possible to recover the exact relativistic result by using 
a A propagator independent of k, even if it reaches the pole (in absence, of course, of the A 
width) inside the allowed region. In fact the present singularity can only be smoothed, even 
in the presence of the A width, after an integration over k is performed, as in the relativistic 
model. 

In general we conclude that, if a dynamical propagator has to be used, like in the semi- 
relativistic approaches discussed here, a smoothed form not hitting the pole is needed in 
order not to deviate too much from the exact result for the response function. Among the 
prescriptions analyzed in this section, the static form of the A propagator produces, in spite 
of its simplicity, the best agreement with the RFG results in the case of one-particle knockout 
sector, once relativistic kinematics and semi-relativistic corrections of the currents are used. 
Therefore, a procedure to "dynamize" the A propagator as the ones presented here, and 
its use within the context of a non-relativistic or semi-relativistic distorted wave analysis of 
quasielastic (e, e') and (e, e'p) reactions appears not to give better results than the static 
approximation. Obviously, this conclusion does not affect the validity of the approximations 
discussed in [28] for the two-particle emission channel, where the kinematical conditions are 
completely different. 

In the next section we apply the present semi-relativistic model of MEC to the exclusive 
(e, e'p) response functions of nuclei. 
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Figure 5: Contribution of the A to the transverse response function using different prescrip- 
tions for the A propagator discussed in the text. 

3 Exclusive (e, e'p) observables 

In this section we compute the quasielastic (e, e'p) response functions for intermediate to 
high momentum transfer values within the context of the SR approach introduced in the 
previous section. As shown in the case of quasielastic inclusive (e, e') responses [30, 31, 
37, 38], accounting for relativistic effects requires at least to treat properly the kinematics 
and relativistic factors in the currents. Referring to the A current, we make use of the 
static limit approach for the A propagator, as this gives rise to the best agreement with the 
fully-relativistic calculation in the inclusive channel, as shown in the previous Section. As 
a complete relativistic analysis of all the MEC (P, C and A) currents in (e, e'p) processes is 
still lacking, the use of "dynamized" A propagators within existing NR or SR descriptions of 
the reaction mechanism does not appear to be well founded and, moreover, these may lead 
to large discrepancies with the exact calculation. 

The general formalism for (e, e'p) reactions has been presented in detail in refs. [2, 3, 4, 45] 
and we refer to them for specifics on the kinematics. Assuming plane waves for the incoming 
and outgoing electron (treated in the extreme relativistic limit) and parity conservation, the 
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exclusive cross section can be written in the form 
d 5 a 



/( = ^<r M (^W L + v T W T + v TL W TL cos 0' + v TT W TT cos 20') , (20) 

where e' and f2' e are the energy and solid angle corresponding to the scattered electron and 
Q p i = (0',(f)') is the solid angle for the ejected proton with four-momentum = (£", p'). 
In eq. (20) K = 2p'mN/(2nh) 3 and a M is the Mott cross section. Finally, v a are the electron 
kinematical factors given in refs. [32, 45]. The labels L and T refer to the longitudinal and 
transverse projections of the current matrix elements with respect to the virtual photon 
direction, respectively. 

The hadronic content of the problem enters via the response functions W a , which are 
obtained by taking the appropriate components of the hadronic tensor 

= 1 £ (p'm s ,^\rm< A) y (p'm.,*?-V\r(Q)\*? ) ) , (21) 

m s M a 

where a sum over undetected final polarization states is performed. In eq. (21), J^{Q) is 
the nuclear current operator and we assume the initial state \&q) to correspond to a spin- 
zero nuclear target in its ground state with energy Eq A \ The final state \p'm s , $^ 4 ~ 1 - ) ) 
is assumed to behave asymptotically as a knockout nucleon with momentum p' and spin 
quantum number m s and a residual nucleus left in a bound state, i.e., I*!?^ 4-1 )) = \ J a , M a ), 
with energy E^ 1 ^. 

In the present distorted wave analysis of (e, e'p) reactions, the outgoing nucleon state is 
described by a wave function solution of the Schrodinger equation with an optical potential 
V opt fitted to elastic nucleon-nucleus scattering data. The matrix elements of the hadronic 
tensor are computed by performing a multipole expansion of both the distorted nucleon 
wave and the current operators in terms of the usual Coulomb Cj, electric Ej and magnetic 
Mj multipoles (see refs. [10, 32, 33] for details on the model, and [17, 24, 32] for explicit 
expressions of the multipole matrix elements of the currents). 

The electromagnetic form factors used in the contact and pionic MEC are given in 
Ref. [24], while for the A we use the same as the electromagnetic proton form factor. 
In this section we include monopole form factors in the MEC F nNN (K) = F wN a(K) = 
(A 2 - m 2 )/(A 2 - K 2 ) with A = 1300 MeV. For low q these form factors are small [24], but 
for high values of q (close to 1 GeV, as here) they can reduce the MEC contribution of about 
25%. 

Note that the A-current contributions can also be sensitive to the A-nucleus potential 
in the medium (see Ref. [21]). Indeed in- medium modifications of the A and ir propagators 
have been studied in inclusive (e, e') scattering in the region of the A peak and for pion 
electro-production [46], but not, to our knowledge, for the MEC in the QEP. A rough 
estimate of these effects can be performed by adding a constant A-nucleus optical potential 
Va = —30 — AOi MeV [47] in the denominator of the static A propagator: since Va turns 
out to affect our results very little, we have accordingly neglected it. 

To illustrate the role of the relativistic corrections we display in fig. 6 the four unpolarized 
exclusive response functions for proton knock-out from the lpi/2-shell in 16 O leading to 
the residual nucleus 15 N. The kinematics selected corresponds to (q, c<j)-constant kinematics 
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q = 460 MeV/c, ui = 100 MeV q = 995 MeV/c, cj = 426 MeV 




p [MeV/c] p [MeV/c; 



Figure 6: Exclusive response functions for proton knock-out from the lpi/2-shell in 16 0. Two 
(q, w)-constant kinematical situations are selected. The traditional non-relativistic results 
(NR) are compared with the responses obtained including relativistic kinematics (RK) and 
using the semi-relativistic (SR) form of the currents. All the curves include MEC. 
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(sometimes also referred to as quasi-perpendicular kinematics) and the values q = 460 MeV/c 
and q = 995 MeV/c have been chosen. In each case the value selected of the transferred 
energy uj almost corresponds to the quasi-elastic peak value, i.e., uj = 100 MeV and 439 MeV, 
respectively. FSI are taken into account for q = 460 MeV/c through the optical potential 
of Comfort and Karp [48], which is appropriate for proton kinetic energies below 183 MeV. 
In the case of q = 995 MeV/c, as the proton kinetic energy is of the order of ~ 430 MeV, 
we use instead a Schrodinger-equivalent form of the relativistic global optical potential of 
ref. [49]. Thus the non- relativistic wave functions correspond to the upper components of the 
relativistic ones, containing the Darwin term. Results obtained within this approach were 
compared with a fully relativistic calculation in the impulse approximation [6], i.e., without 
including the effects of MEC. 

In each panel of fig. 6 we compare the results corresponding to the traditional non- 
relativistic model with non- relativistic kinematics (denoted by NR), including relativistic 
kinematics (RK), and finally, the semi-relativistic approach discussed in the previous section 
(SR). All the calculations include MEC. From these results it emerges that, particularly for 
high q, the relativistic kinematics plays a crucial role in describing properly the form of the 
momentum distribution. Note that the allowable missing momentum values are determined 
by the relation \p' — q\ < p < p' + q, with p' fixed by energy conservation. Assuming 
non-relativistic kinematics the value of p' is given by 



p' = ^/2m N (e h + Lu), (22) 
while for relativistic kinematics it results 



p' = J (m N + e h + uj) 2 - m 2 N = J 2m N (e h +u)[l + ^ — - ) . (23) 



2m 



N 



In eqs. (22,23), the energy represents the (negative) energy of the bound nucleon. Hence in 
the case of RK we solve the Schrodinger equation with equivalent energy (e^ + uj) (l + f^~ ) 
for the ejected proton instead of the NR energy eh + uj. 

Once relativistic kinematics is selected, the use of SR currents produces an enhancement 
of the L and TL responses and a reduction of the T and TT ones, whose magnitude increases 
with q. These effects are connected with the factors -7= > 1 and ~ -p-— that enter in 

^ y/r k 1+r 

the SR expressions of the longitudinal and transverse currents, respectively [31, 37]. It 
is important to point out the particularly large relativistic enhancement observed for the 
interference TL response. This is connected to the spin-orbit term [32] contained in the SR 
charge density but neglected within the NR approach. As shown in [32], the interference 
between the spin-orbit term and the magnetization current gives rise to a contribution in 
W TL which is of the same order of magnitude as the interference between the charge density 
and convection current. Thus the presence of the spin-orbit term in the density is essential 
to describe properly the response W TL , even for moderate q- values. 

Analogous results hold for a proton knockout from the p 3 / 2 -shell in 16 O: the role of 
relativity in each response is similar to the one observed for the pi/2 case and, again, the 
most sensitive response to relativistic effects is W TL . However, the relativistic enhancement 
for the p 3 / 2 -orbit is somewhat smaller than the one presented for pi/ 2 , particularly for large 
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Figure 7: Separate contribution of each one of the MEC to the exclusive response functions. 
SR results are shown for knock-out from the lpi/2 shell of 16 for the same kinematics as 
in fig. 6, including the one-body current only (OB), and in addition the contact (C), pionic 
(P) and A currents in the calculation. 
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Figure 8: The same as fig. 7 for knock-out from the shell of 16 0. 

The separate contributions of the MEC is presented in figs. 7 and 8 for p 1 / 2 and ps/2- 
shells, respectively. Kinematics is as in the previous figure and we do not show results for the 
pure longitudinal response as it is not affected by MEC within the present approaches. We 
compare the OB responses (dotted line) with the results obtained when including the contact 
(C) current (dashed line), the contact and pion-in-flight (C+P) currents (dot-dashed line), 
and the contact, pionic and A (C+P+A) currents (solid line). Let us analyze each response 
separately. In the case of W T , the contact current produces an increase which is partially 
cancelled by the tiny reduction introduced by the pionic current. Note that the role of the C 
and P currents is negligible for large q. The A current gives rise to an additional reduction 
whose relative magnitude is rather similar for both shells. The total effect of C+P+A MEC 
is a reduction of the T response which is slightly more important for high q. 

Although a great caution should be taken in extending the conclusions drawn from the 
analysis of exclusive reactions to inclusive ones or viceversa, it is illustrative to discuss the 
results relative to quasielastic (e, e'p) reactions in connection with the quasielastic (e, e') 
ones. Note that, apart from the potentials usually considered in both types of reactions, the 
inclusive responses are calculated by integrating the exclusive ones over the ejected nucleon 
variables, and summing over all occupied hole states, i.e., including also the contribution 
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given by the lsi/2 shell. In refs. [31, 37, 38] we evaluated the role of MEC for the inclusive 
T response in a relativistic Fermi gas model. The contact (C) contribution was found to 
be larger than the pion- in- flight term (P), a dominance increasing with q. The total (C+P) 
contribution presents an oscillatory behaviour with respect to the transferred energy u>, 
having a node close to the QEP value, particularly for high q. For q — 1 GeV/c, this node 
is reached for u ~ 435 MeV, which is very close to the value oo = 439 MeV selected in the 
(e, e'p) calculations (right hand panels in figs. 7 and 8). On the contrary, in the case q = 500 
MeV/c, the contribution of the C+P currents in the inclusive T response presents a node 
for uo ~ 150 MeV and is positive and non-negligible for the value oo = 100 MeV selected here 
in the case of (e, e'p) process. A similar discussion can be also applied to the A results. In 
ref. [38] (see also previous section) the A contribution to the inclusive T response has been 
shown to be always negative up to q ~ 2 GeV/c and larger than the one provided by the C 
and P currents together. Hence the general trend of the MEC contribution to the inclusive 
channel agrees with the effects shown for the exclusive responses (figs. 7 and 8). 

Finally, comparing the results in figs. 7 and 8 we conclude that the effects of MEC for 
the exclusive T-response are similar for both pi/ 2 and p 3 / 2 shells. A different behaviour is 
found in the results of ref. [25] where the role of MEC, for q = 460 MeV/, is shown to be 
larger in the case of the p 3 / 2 shell. In addition, the role introduced by the pionic current 
(P) in [25] is said to be negligible, whereas in our case its contribution, though a little bit 
smaller than the contact term, is clearly visible. Our results are also in disagreement with 
the calculation of [21] performed for q = 1 GeV, where the total MEC effect is found to be 
small for the p\/ 2 and large for the ps/ 2 in the T response. This discrepancy is linked to the 
FW method used to implement relativity in [21] which is not expected to be applicable for 
g=l GeV [37]. 

Next we focus on the interference TL response. From figs. 7 and 8 we observe that for 
q = 460 MeV/c the MEC effect is larger in the case of the pi/2-shell. The main contribution 
comes from the contact current, which produces an enhancement in W TL of the order of 
~ 20% for pi/2 and ~ 15% for ps/ 2 , while the role of the pion-in-flight is smaller, reducing 
the response, and in particular, the A gives rise to an almost negligible contribution (slightly 
positive for p X j 2 and negative for p 3 / 2 ). These results disagree with the findings of ref. [25], 
where the effect of the A in W TL is substantially large and negative for the p^-shell, so 
that it cancels out the contact term, yielding a negligible global MEC contribution. On the 
contrary, in the same reference the A contribution for the ^3/2-orbit is found to be large and 
positive, hence the net effect of MEC is an important enhancement in the TL response. As 
in the case of the pure T response, the pionic contribution found in [25] is negligible, which 
again is not in accord with our results. 

Figs. 7 and 8 clearly show that the importance of the MEC on the TL response decreases 
as q goes to higher values. Whereas the contact current enhances the TL response by the 
same magnitude for the two p-shells, the A term is negligible for the p\/ 2 and tends to cancel 
the contact contribution for p s / 2 . For both shells, the pionic current does not alter W TL . 

As far as W TT is concerned, note that this response is much smaller, its contribution being 
of order {kp/m^) 2 (see ref. [32]); hence terms of second order in p/m^, usually neglected in 
the expansion of the current operators, may provide a crucial role in this response. However, 
since W TT is also particularly sensitive to the details of the model, it can be used as a test to 
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Figure 9: A^L-asymmetry for proton knock-out from the lpi/2 and lp^/2 shells in 16 0, 
computed within the SR model including the successive contribution of each one of the MEC. 
Experimental data are from [8]. Left and right panels have been obtained to including or 
not the spin orbit part in the optical potential. 

compare different theoretical models. Note that this response is opposite in sign for p 1 / 2 and 
P3/2- Concerning the role of MEC, it is found to be of the same relative order of magnitude 
as in the other responses, except for the pionic current which has an important effect in this 
case. 

As a further application of our model we present in fig. 9 the TL asymmetry (A TL ) for py 2 
(top panels) and p 3 / 2 (bottom panels). The kinematics has been selected to correspond to the 
experimental data [8] . Atl is obtained from the difference of cross sections measured at 0' = 
0° and 0' = 180° divided by the sum, hence this observable is particularly interesting because 
it does not depend on the spectroscopic factors. A detailed study on the TL asymmetry has 
been carried out in refs. [6, 7] within the relativistic distorted-wave impulse approximation 
(RDWIA), namely a fully relativistic calculation without including MEC. In particular, a 
comparison between A TL evaluated in RDWIA and using the SR approach for the one-body 
current and neglecting MEC, was presented in [6]. There it was proved the large contribution 
given by the spin-orbit correction to the charge density, and, more importantly, the crucial 
role played by the dynamical enhancement of the lower components of bound Dirac spinors 
in the description of interference TL observables (see also refs. [50, 51]). Recent data on 
polarization observables nicely agree with RDWIA analysis [52, 53]. 
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However, a direct calculation of the role of MEC in Atl within a relativistic approach 
has never been performed, as far as we know. Thus, in fig. 9 we show the results obtained 
within the impulse approximation, i.e., without MEC (dotted line), and the contributions 
introduced by the various two-body currents: C (dashed line), C+P (dot-dashed line) and 
total MEC, namely C+P+A (solid line). In order to investigate also the dependence of our 
results on the FSI, in fig. 9 we show two sets of calculations for each p shell: including the 
spin-orbit part of the optical potential Vi s [49] (left panels) and without it (right panels). For 
comparison we also present the experimental data [8] . From inspection of fig. 9 we conclude 
that the effect of MEC is very small for low missing momentum values and it starts to be 
important for p > 300 MeV/c. Note however that in this region the dynamical enhancement 
of lower components, obviously not considered within the present SR approach, starts also 
to play a crucial role and hence it should be considered before a detailed comparison with 
data can be accomplished. Note also the large discrepancy between the SR calculation and 
the data for low p in the case of p±/ 2 with the full potential. This issue was already present 
in [6], where the RDWIA calculation differs from the SR one comparing better with data [8]. 
This problem within the SR approach appears to be connected with the spin-orbit term 
introduced by the equivalent Schrodinger form of the optical potential, as can be seen in the 
right panels of fig. 9, that do not include Vi s in the FSI. Our results also show that the MEC 
effects for high missing momenta strongly depend upon the FSI, since they are substantially 
reduced when only the central part of the optical potential is included (right panels). 

Summarizing, from theoretical results in fig. 9 and ref. [6], and their comparison to 
experimental data [8], we may conclude the following: i) a fully relativistic calculation 
within the impulse approximation (RDWIA), i.e., including the effects introduced by the 
dynamical enhancement of the lower components in the Dirac spinors, appears to be essential 
to reproduce the data; ii) the effects introduced by MEC for high missing momentum values 
seem to be also very important, and highly dependent on the FSI. Hence an appropriate 
relativistic analysis of these two-body currents may be also essential in order to improve the 
description of the experimental data at high missing momentum. 

4 Conclusions 

In this paper we have presented a semi-relativistic model of inclusive and exclusive elec- 
tron scattering from nuclei in the one-nucleon emission channel, including one- and two- 
body currents. These currents differ from the usual non-relativistic ones by multiplicative 
(q, u;)-depending factors obtained by an expansion in powers of the missing momentum. An 
essential ingredient of the SR consists in using, in addition to the SR currents, relativistic 
kinematics to relate the energy and momentum of the ejected nucleon. 

This model has been already tested in quasielastic inclusive (e, e') processes in refs. [31, 
37, 38]. In this paper we have first focussed on the new SR-A current and compared its 
contribution to the inclusive transverse response with a fully relativistic calculation in RFG 
for moderate and high g-values. This has allowed us to test the reliability of different 
prescriptions introduced to account for "dynamical" aspects of the A propagator. We have 
found that the best agreement with the fully relativistic calculation corresponds to the 
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static limit approximation. Thus, we conclude that the use of any of these "dynamized" A 
propagators within the one-particle emission channel is not justified and, moreover, it may 
produce very large discrepancies with the exact result. 

Next we have implemented the SR currents into a DWIA model of the quasielastic (e, e'p) 
reaction, computing the separate response functions and the TL asymmetry, using the static 
limit for the A propagator. After analyzing the role of the different relativistic corrections 
embedded in our calculation and its dependence upon the momentum transfer, we have 
studied the effect of MEC on the different observables. In particular, we have compared the 
results for proton knock-out from the pi/ 2 and p 3 / 2 shells in 16 0, choosing quasi-perpendicular 
kinematics, typical of the experiments. In the case of the T response, MEC effects are shown 
to be equally important for the two shells, giving rise to a net reduction of the response. For 
W TL , the role of MEC is to enhance the response, an effect that is substantially larger for the 
P1/2 shell, and reduces considerably for high q- values. In general, we get sizable differences 
with previous calculations in the literature [21, 25]. 

Finally, our model applied to the A?l asymmetry shows that this observable is very 
sensitive to the MEC and to FSI in the region of high missing momentum. Since it has been 
proved in [6] that the interference TL observables, particularly Atl, are also crucially affected 
by other relativistic ingredients, such as the dynamical enhancement of lower components, 
not included in this work, it would be very interesting to evaluate MEC effects within the 
scheme of the fully relativistic calculation of refs. [4, 6, 7], and contrast their predictions 
with the ones obtained with the present SR model. 
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